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The extraction of the finite temperature heavy quark potential from lattice QCD relies on a 
spectral analysis of the Wilson loop. General arguments tell us that the lowest lying spectral peak 
encodes, through its position and shape, the real and imaginary part of this complex potential. Here 
we benchmark this extraction strategy using leading order hard-thermal loop (HTL) calculations. 
I.e. we analytically calculate the Wilson loop and determine the corresponding spectrum. By fitting 
its lowest lying peak we obtain the real- and imaginary part and confirm that the knowledge of the 
lowest peak alone is sufficient for obtaining the potential. Access to the full spectrum allows an 
investigation of spectral features that do not contribute to the potential but can pose a challenge 
to numerical attempts of an analytic continuation from imaginary time data. Differences in these 
contributions between the Wilson loop and gauge fixed Wilson line correlators are discussed. To 
better understand the difficulties in a numerical extraction we deploy the Maximum Entropy method 
with extended search space to HTL correlators in Euclidean time and observe how well the known 
spectral function and values for the real and imaginary part are reproduced. Possible venues for 
improvement of the extraction strategy are discussed. 



I. INTRODUCTION 

Twenty seven years ago Matsui and Satz [1^ proposed 
the melting of J/^, i.e. the ground state of the cc vec- 
tor channel, as signal for the deconfinement transition 
in heavy-ion collisions. The recent success of relativis- 
tic heavy-ion experiments [2H5] in observing the relative 
suppression of charmonium and bottomonium serves as 
further motivation to develop a first principle description 
of the phenomena. 

In the framework of effective field theories, heavy 
quarks can be described by non-relativistic quantum 
chromodynamics (NRQCD) obtained form QCD by in- 
tegrating out the hard energy scale, given by the rest 
mass of the heavy quarks. To describe the bound state 
of two quarks, one can further integrate out the typi- 
cal momentum exchange between the bound quarks (see 
[6] and references therein), which leads to potential non- 
relativistic QCD (pNRQCD). In this effective field theory 
the bound state is described by a two point function sat- 
isfying a Schrodinger equation. 

At zero temperature, the potential between a heavy 
quark and anti-quark is defined from the late time be- 
havior of a Wilson loop and can be directly calculated 
in Euclidean-time lattice simulations or in perturbation 
theory. At small distances, where perturbation theory 
converges, both results agree J7J. 

At high temperature, above the QCD phase transition, 
one might first expect that the problem becomes simpler 
as the potential is not confining anymore. Actually, this 
is not the case since even a proper definition of the poten- 
tial becomes non-trivial. In fact, the presence of a heat 
bath is most conveniently incorporated in a Euclidean 
time framework with finite temporal extend. There the 
Wilson loop depends on imaginary time and needs to be 
analytically continued to real time. Only from the large 



real-time, i.e. i — >■ oo behavior, the finite temperature 
potential can be extracted and happens to be complex 
[H E] ■ Its imaginary part can be interpreted as Landau 
damping [10] and describes the decaying correlation of 
the QQ system with its initial state due to scatterings in 
the plasma. 

Along the lines presented in [5], one can compute the 
potential in finite temperature perturbation theory. This 
is a demanding task, as resummations need to be carried 
out in order to cure infrared divergences. To this day the 
full result is known only to leading order, whereas a short 
distance expansion has been calculated to higher order 
[TTl [T^ . Even if higher orders were available, observing 
the deconfining transition will remain beyond the reach 
of perturbation theory. 

In Ref. [13], a method was proposed to compute the 
heavy quark potential non-pcrturbatively from lattice 
QCD simulations. Starting from the measurement of the 
Euclidean Wilson loop on the lattice, its spectral func- 
tion is reconstructed via the maximum entropy method 
(MEM). The definition of the potential is based on the 
peak structure of the Wilson loop spectrum. 

Previous numerical evaluations however lead to unex- 
pected results: both the real and imaginary part appear 
to grow linearly at distances where other quantities, such 
as the free energies already show significant screening ef- 
fects. This behavior persisted even at temperatures much 
larger than the QCD phase transition, where on general 
grounds, one would expect that the confining potential 
disappears because of Debye screening [14] . 

This problem was solved recently |T5] by carefully dis- 
entangling the different timescales in the problem. Tak- 
ing into account the remnants of early-time non-potential 
physics, the lowest lying spectral peak was found to devi- 
ate from a naive Lorentzian shape through skewing. Ex- 
tracted values for real- and imaginary part based on this 



functional form result in a potential that is compatible 
with Debye screening. 

In this paper our aim is twofold: at first we wish to as- 
certain whether fitting of the lowest lying spectral peak 
indeed suffices to determine the static heavy quark po- 
tential, given the spectral function of the Wilson loop 
or even the gauge fixed Wilson line correlators. Subse- 
quently it is our goal to better understand the challenges 
facing a numerical determination of the spectral function 
by Bayesian analytic continuation. Since in the pertur- 
bative approach both Euclidean correlator and spectrum 
are known, the outcome of the numerical reconstruction 
can be readily compared. 

In section [ill we review the basics of the method of 
Ref. ,13j and its improvement introduced in [T5 , which 
form the basis of the extraction of the potential from lat- 
tice simulations. From calculations of the real-time Wil- 
son loop as well as gauge fixed Wilson line correlators in 
section [III] we determine and investigate the correspond- 
ing spectral functions in section IV While in section[V[we 
apply the peak fitting procedure of ^liji to the HTL spec- 
tra, section \VJ\ scrutinizes how well these spectra can be 
obtained with the maximum entropy method from the 
HTL Euclidean correlators. Our conclusion in section 



VII discusses the limitations of the method and points 
toward further possible improvements. 



II. HEAVY QUARK POTENTIAL FROM 
EUCLIDEAN CORRELATORS 



The description of the interactions between a pair of 
heavy quarks and antiquarks at finite temperature in 
terms of a quantum mechanical potential V(r) requires 
the relevant physics to be well separated from the energy 
scale of pair creation. In particular 
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needs to be fulfillecQ which is satisfied exactly in the 
static limit (toq — ^ oo). In that case, the propagation 
amplitude of an infinitely heavy quark pair can be de- 
scribed by a rectangular temporal Wilson loop Wa{t,r) 
where t, r are its temporal and spatial extend. This real- 
time quantity is defined as the closed contour integral 
over the matrix valued gauge field v4^(a:) = ^^(x)T° 
along the path of the heavy quarks 



Wait,r) 






)]]• (2) 



If the scale hierarchy holds, it is permissible to substitute 
the field theoretical interactions by an instantaneous po- 



tential, so that W\j{t,r) obeys a Schrodinger type equa- 
tion 



idtWait,r) = ^t,r)Wa{t,r) 



(3) 



At late times, on expect the function $(t, r) to become 
time independent, so that we may define the heavy quark 
static potential as 



V{r) = lim $(/ 



t). 



(4) 



Due to the complex weighting factor in Feynman's path 
integral, we cannot calculate the real-time Wilson loop 
using lattice QCD Monte Carlo simulations. Instead we 
have to rely on an analytic continuation of Euclidean 
time quantities that are accessible by these numerical 
methods. In order to connect the heavy-quark potential 
V{r) and the Euclidean Wilson Loop one introduces a 
spectral representation of the real-time quantity. 



(5) 



Wair,t)^ / du e-'^' pair,u;), 



where the time dependence now resides entirely in the in- 
tegral kernel. Note that the function pai^T^) is not just 
a Fourier transform but can be shown to be a positive 
definite spectral function fTTrj After analytic continua- 
tion t = —IT one observes that only the integral kernel 
has changed, whereas the spectral function remains the 
same 



Wa{r,T)^ / rfc^e-">n(r,w). 



(6) 



Using the Maximum Entropy Method (MEM) , a form of 
Bayesian inference, it is in principle possible, albeit chal- 
lenging, to invert Eq.(|6]) and thus to extract the spectral 
function from Wn (r, r) . Note that the model indepen- 
dent method of refs. [52H23] is not directly applicable 
as the Wilson loop is not periodic. However a similar 
method could probably be developed from the general 
results of refs. [251 [251. Once we are in possession of the 
spectral function /?□ we can insert Eq.((5]) into Eq.Q, 
which yields [17] 



V{r) — lim 



duj LU e-*"^* p{r, u) 
duj e-'"^^ p{r,u) 



(7) 



Direct application of this formula in the case of a numer- 
ically reconstructed spectral function is very difficult. It 
is however possible to determine those structures in the 
spectral function, which dominate the integral in the in- 
finite time limit. 



^ See for instance |16| for the discussion of the different limiting 
cases and their physics 



^ It is important to distinguish this r-dependent Wilson loop spec- 
tral function from the quarkonium spectral function |18H21j rep- 
resenting the physical quarkonium spectrum. 



If we suppose that the time independent potential de- 
scription holds for all times t i.e. <&(t,r) = V{r) in equa- 
tion (pi, an intuitive connection between spectral features 
and the static potential can be established. In this case 
equation ([3]) can be solved and the spectral function turns 
out to be a simple Brcit Wigner peak 



piuj,r) 



lin[V]{r) 



lin[V]{r)^ + {Re[V]{r) ~ Lj) 



(8) 



characterized by its peak position ^0(7") — Ke[V]{r) and 
width ro(r) ^lm[V]{r). 

In general the function ^{t,r) however is time depen- 
dent at early times and one expects that a wealth of 
structures, different from the simple Lorentzian exam- 
ple, exists in the spectrum of the Wilson loop at finite 
temperature. Note that if the potential description is ul- 



timately applicable, the function <I>(r, t) will become time 
independent at late times and therefore a corresponding 
well defined lowest peak must exist. This part of the 
spectrum encodes all the relevant information on the po- 
tential and it alone needs to be reconstructed from the 
Euclidean correlator. 

In Ref. [13! it was assumed that the lowest peak is 
solely described by the late time behavior of the poten- 
tial and is not affected by the time dependence of the 
potential at short times. It was shown in Ref. [15] that 
this is actually not the case. The short time dynam- 
ics (non-potential terms, bound state formation) doesn't 
just create additional structures at high frequency but 
also significantly modifies the shape of the low frequency 
peak. The most general form of this low peak, derived in 
Ref. [15], can be written as 



PD{r,uj) = -e 

TT 



M 



,) |Im[y](r)|cos[Re[aoo](r)] - (Re[F](r) ~ w)sin[Re[aoo](r)] 



Im[F](r)2 + (Re[y](r)-w)2 
+co(r) + c,{r)tQQiRe[V]ir) - oj) + C2{r)tl^{Rc[V]{r) - uf 



(9) 



Note that this result can also be obtained from pNRQCD 
where Re[cr 00] arises from the phase of the singlet nor- 
malization factors Zs (r) [6] 

In order to calculate the potential V{r) from Euclidean 
correlators we thus need to carry out the following steps: 

1. Calculate the Wilson loop Wui^TT) at several sep- 
aration distances r for all possible values along the 
imaginary time axis r g [0,/3]. 

2. Use Bayesian inference to extract the most proba- 
ble spectrum pa{r, 1^) for each value of r. 

3. Use Eq. (It]) to determine the potential V{r) 

(a) by direct Fourier transform of the full pn (^j i^) j 
which is usually impractical due to the uncer- 
tainties introduced by the MEM OR 

(b) by fitting the lowest lying peak with the func- 
tional form (|9J) and analytically carrying out 
the Fourier transform in Eq. ^ 

In the following section [llj we prepare a testing ground 
for this extraction strategy based on analytic calculations 
of the real-time and Euclidean Wilson loop in the HTL re- 
summed perturbative approach. Since the analytic con- 
tinuation can be performed explicitly in HTL, item three 
of the above list can be tested independently from ques- 
tions arising from possible inadequacies of the maximum 
entropy method. The availability of both the spectrum 
and Euclidean data points on the other hand furthermore 
allows us to check the degree of success of the MEM itself 
in the form of a realistic mock data analysis. 



III. CORRELATORS FROM HTL RESUMMED 
PERTURBATION THEORY 

A. Wilson loop 

In perturbation theory, the Wilson loop is calculated 
as an expansion in the coupling: 

Wa{r,T) = W^\T,r)+g^W^\r,r)+0{g^), (10) 

starting form W^ =1. The first non-trivial term [W^ ) 
contains only a one gluon exchange and is not enough to 
describe the correct physics for large Euclidean time r. 
To improve this situation, we resort to the usual 'expo- 
nential' resummation [TU], noticing that 

\og{Wa{T,r)) = g^W^\T,r)+0{g^). (11) 

Thus a better approximation for VFn {^i f) is 

Wu{T,r)^eM9^wg\r,r))+0{g^), (12) 

as it resums all 'ladder diagrams' and contains the correct 
leading order (g^) large r behavior. 



1. Leading order term 

(2) 
We now turn to the calculation of W^ (t, r), for which 

we set the r direction along the third spatial axis. In hard 

thermal loop (HTL) resummed perturbation theory, all 



(2) 

diagrams contributing to W^ (r, r) have one HTL gluon 
running between the Hues of the Wilson loop [8J : 



W^'^ = CfT 



-iqoT-i 



(27r)3 2 

T2Aoo(0,g)+^(2-e'««- 

2Ao3(go,g) , A33((7o,9) , Aoolgo,^) 



(13) 
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The gluon HTL propagator, written in Euclidean space 
{Q^ — Qi + Qq) a-iid covariant gauge reads: 



A^.(Q)=rt e^'^^^-y) 



Q^+TIt{Q) 



PLiQ) 



iti^qv 



Q^ + iidQ) '(Q2)2 



(14) 



while the HTL self-energies TiE,T are given in Appendix 
[X|and the projectors take the form: 



PmiQ) - PliQ) = PfoiQ) = 0, PfAQ) = <5, 



PLiQ) = s,. - 



Q' 



-PLiQ)- 



qiQj 



(15) 



Following Rcf. [S], we rewrite the HTL self-energies as 
spectral functions, 



1 

g^+q2 + HL,T(9o,q) 



"-M^Ei^, (10) 



oo 



TT q" 



na 



so that we can perform the sum over go analytically: 



W^\T,r) = CF 



(27r)3 

T 



+ e-"'^'' - 2 



(17) 



q2+ni(0,q) 



PL(g°,q) 



dgO 



nB(g°)M^,'?°) 



(9") 



o^2 



+PT(<7°,q)f;^-^ 

\y3 '^ 

where we abbreviated the r dependence of the second 
term through the function 



MT,gO) 



1 



Ji" 



Af>-r)q° 



, grq _ ^(p-r)q (18) 



We can write the spatial vector q in spherical coordinates 
(9 = IqIi^:*^) a-nd 93 = qcoaO. In an isotropic plasma, 
the HTL spectral functions and self-energies depend on 
q, q° only. Integrating over (j) is trivial and the integral 
over c = cos 9 involves 



^'^^ + ^"^"^^dc = 2f^^^-l). (19) 
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dc = 2 (1 — cos(<7r) — rq'&i{qr)) 



where Si is the sin integral function. Performing the an- 
gular integrals and using H£;(0, q) — rn^ gives 

"^ ^ sva{qr) 



W'^>{T,r)^CF 



dq° 



dq_ g' 
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(20) 



2- 



sin(gr) 



gr 



3(gr) -grSi(gr) j pT(g°,g)|- 



The first line of equation (21) is linear in r, whereas 
the next lines are proportional to /i(t, g°) and therefore 
symmetric around r = 13/2. We will consider these terms 
separately in the following: 

W^\T,r) = W^£^(r,r) + <l(r,r). (21) 



2. Part linear in t 

The part linear in r is formally divergent. Using di- 
mensional regularization, the result can be read off from 
Ref. (H|; the first line of equation (21) hence gives: 



wl^iT,r)^Cf 



dq 

2^ 



sin(gr) 



-1 



qr 



tCf 

47r 



mD 



(22) 



In the limit t ^^ it ^ ioo this part yields the real part 
of the potential: 



d 



Re[V]{r)^g\]imi-m^{zt,r) 



rm, 



t->-oo dt 

in 



(23) 



mD 



Note that the result is finite (for r ^ 0) and the di- 
vergence at r = reflects the behavior of the Coulomb 
potential. 

On the lattice, this term behaves differentljjj Roughly 
speaking, the integral is truncated by the lattice cutoff 
g < A and thus finite. In this case it is easy to see that 
it vanishes at r = 0, which is expected as a Wilson loop 
without area is equal to unity. For r > 0, it decreases 
quickly and formally goes to — oo in the limit of an infinite 
cutoff. 

This behavior cannot be canceled by the other terms 



in equation (21) as they have a different r dependence. 



It should also not be removed as it encodes the Coulomb 



^ The difference with dimensional regularization can be traced 
back to an infinite constant that is removed in the dimensional 
regularization procedure. 



part of the potential that we want to obtain. To make 
a connection to the lattice, we therefore introduce a UV 
cut-off, mimicking the finite lattice spacing. In this case, 
performing the integral over the momentum q from zero 
to A in equation (22) gives: 



the pole contribution of the electric spectral function as: 



XTr(2),A. ^ ^ T 



—A + Tod tan 



A 
Tod 



cosh(?TiD7')(Si(r(iTOD — A)) — Si(rA + irTO,D)) 

(tt — iCi(rA — irrriD)) + iGi{irmY) + Ar)) sinh(7TtD'') 



where Si, Ci are the sin and cos integral function. From 
the UV regularized version of the correlator we get the 
following potential. 



d 



^(2)>A/,, 



Re[^^](r) = g^ lim ^^M^//J'^^(^i, r), (24) 

which is plotted in Fig. [7] together with the continuum 
(A — > oo) potential. 



3. Symmetric part 



We calculate here the symmetric part of the correlator 



.(2) 



Wsym{T, r) containing the lines 2-4 of equation (21 1. The 



functions PL,T{,<f , l) receive a contribution from the cuts 
of Ul,t if (? > \q^\- For the opposite case |g°| > q they 
vanish except for a (^-function contribution coming from 
the pole of n^.T- In the following we calculate the con- 
tribution from the cuts and poles of the transverse and 
longitudinal self-energy separately. 



" sym "cut ' " pole,L ' " pole.T 



(25) 



As before we introduce a cutoff on the momentum to 
mimic the effects of the lattice regularization. 

a. Cut contributions Using the symmetry go ^ 
—qo, the cuts contribute to the Euclidean Wilson loop 
as 



W^lliT,r) = CF 



"^ dq /"^ dq° 
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cos(gr) ~qrSi{qr)j Pt(<?",9)|, 



where the integrals should be performed numerically and 
the functions pl.t are given in Appendix [A| Note that 
in Eq. (261, the limit A -^ cxo is well defined. 

b. Pole contribution form the longitudinal spectral 



function We can write the part of (21| coming from 



.(2) 



KoLLir,r) = C, 



dq 



dq"nB{q'')HT,q") 



'sin(gr)_ 



qr 

A 



= Cf I -^nB{ql)h{T,ql) 

Jo 7^^ I/l(9l)I 

Here q^ j, is the solution of fL,T{Q^) — 0, g° > and the 
remaining integral is performed numerically. The limit 
A — > cx) also exists in this case (see Appendix IbI) . 

c. Pole contribution from the transverse spectral 
function We proceed in a similar way for the transverse 
spectral function. 



KolTir,r)^C, 



dq . ^ H'T, Qt) 



r"B(9T)T77 



\M'i"t)\ 



(28) 



/„ sin((7r) , , „./ N 

X 2 cos(gr) — qrbi(qr) 

Here, the limit A — )> cx) does not exist the integral in equa- 
tion (281 is linearly divergent (see Appendix Ib| . Note 
that such divergences were already observed iii[27l |28| , 
where the Wilson loop of maximal time extend t = (3 
is shown to diverge at next to leading order. The lead- 
ing order divergence found in Eq. [28] has yet a different 
nature and consistently vanishes for r = 0, /3. In dimen- 
sional regularization, it can be shown (see appendix O) 
to match the cusp divergence [351 ISO] : which in this case 

gives "^m- ^ 

Here, we are not interested in trying to renormalize 
the Wilson loop. It is not needed for our purposes as we 
aim at a comparison with lattice results, which are also 
not renormalized. It is however interesting to note that 
these cusp divergences do not contribute to the potential 
and only make the Wilson loop heavily suppressed for 
T 7^ 0, /3, hence harder to measure with high accuracy. 
Removing these divergences in the lattice measurements, 
without affecting the potential would be of great help to 
improve the accuracy of the lattice data. One strategy 
deployed to this end could be the smearing of gluonic 
links [S]- 



4. Imaginary part of the potential 

From the symmetric part, we obtain the imaginary 
part of the potential, 

zlm[y]^(r) = g^ lim i^W^H^iit, r). (29) 

As in the end the infinite time limit will be taken, it is 
sufficient to consider the low frequency part of the if 



integrals, 



MV]\r)^9' lim l-l % r^n^(,0)-| 
^ ^ ^ > ^ t-^^dtJo TT^ Jo TT """-^ \qO 
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_o I sm{qr) 



qr 
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Performing the time derivative, using equation (A3) and 



approximating nB{q^) ~ T/q^ for small q^ as well as the 
identity 

itq° _ Mi-it)q° 

lim -^ = 2TTi5{q'^), (30) 



we get: 



lm{V]'^{r) = - 



47r 






which coincides with the expression obtained in [SUlOj. 



5. Numerical evaluation 

To make close connection to actual lattice data with 
spatial lattice spacing a ~ 0.04fm, we choose to fix the 
cut-off in our HTL calculations to 



A = 



(31) 



which naively corresponds to the largest momentum ac- 
cessible under this finite resolution. Based on a numerical 
evaluation of the remaining integrals in eq. ( 2f|25p8 ), 
we can generate an arbitrary large number of datapoints 
spanning the imaginary time axis, which carry numerical 
errors of the order of the machine precision only. 

Comparing this ideal HTL Euclidean regularized data 
to actual measurements from a Monte Carlo simulation 
in Fig. [l] we find a strong qualitative resemblance. Both 
graphs exhibit three characteristic features, i.e. a sup- 
pression region at small r together with an upward trend 
at r ~ /3, both of which are closely linked to the diver- 
gences observed in |III A3"| The datapoints at intermedi- 
ate r are the ones encoding the potential. They exhibit 
nearly exponential behavior for small separation r, where 
also Im[y] is small but begin to show noticeable curva- 
ture for larger separation distances. 

After calculating the real-time values Wu{it:r) (see 
Fig. ^ using a similar numerical evaluation of the inte- 
grals in ( 21|25p8 ), it is possible to obtain the function 
<i>(t, r). As shown in Fig. [2J we can explicitly observe 
the approach of $(i, r) to a constant value and thus the 
emergence of a simple exponential behavior of the Wilson 
loop. Note that in Fig. [2] we show times t < 40GeV^^ 
where the oscillatory behavior is clearly visible while a 
constant value is actually reached for larger t. We re- 
frain from attaching any physical meaning to the length 
of the swing-in period, as it is dominated by the same 
cusp divergences that lead to the suppression of the Eu- 
clidean Wilson loop data points. 



!□ 



1e+00 
1e-01 
1e-02 
1e-03 
1e-04 
1e-05 
1e-06 
1e-07 
1e-08 
1e-09 
1e-10 



1e+00 



1e-01 



1e-02 



1e-03 



!••-••#, ***••••• ••••• 

••• *••• ***••••• ••••••• 

•**•**•• **••••• •••••« 

Ar=0.066fm *•••.-,- „ , c • * 



10 



15 



20 



25 30 



1e-04 



P=7.0 
Ar=0.039fm 



'••i. 



o 

• • • » • T 

..J.I 

. • . • 







10 



15 



20 



25 30 



FIG. 1. (top) The Euclidean HTL Wilson loop Wa^^(T,r) 
with momentum regularization A = Stt GeVevaluated at T = 
2.33 X 270 MeVin steps of Ar = 0.066fm. (bottom) Quenched 
lattice QCD Wilson loop from a lattice with a = 0.039fm and 
anisotropy ^ = 4 at T = 2.33rc. 



B. Gauge fixed Wilson line correlator 

Cyclic Wilson line correlators (i.e. color singlet 
Polyakov loops) fixed to Coulomb gauge have been exten- 
sively studied on the lattice, both for the determination 
of the zero temperature potential as well as in investiga- 
tions into the free-energy difference between a medium 
with and without inserted heavy quarks (see for instance 
[32l [33] ) . Due to the absence of spatial Wilson lines con- 
necting the temporal links, these quantities offer a signif- 
icantly better signal to noise ratio than the Wilson loop, 
especially if the multilevel algorithm [34j is applied. 

Besides the technical question of whether the removal 
of spatial connectors (or e.g. the application of smearing 
on spatial links) can lead to an improved lattice observ- 
able for the extraction of the potential, it is conceptually 
of interest to understand whether gauge independent in- 
formation, such as the potential can be extracted from a 
gauge dependent quantity such as the Wilson line corre- 




FIG. 2. (top) The HTL real-time Wilson loop Wg'^^{t,r) 
with momentum regularization A = tt GeVevaluated at r = 
IGeV"^ and T = 2.33 x 270 MeV. (bottom) Time evolution of 
the quantity ^{t,r) obtained from W§'^^{t,r) through EqM 



latorqf] 

We proceed with the determination of the Euclidean 



time Wilson line correlator analogously to III A 

W||(r,r) = l+.g2l^lf)(r,r) +©(/). (32) 

Calculating in leading order hard thermal loop (HTL) 
rcsummed perturbation theory we obtain the expression 



|y(2) ^ (j^j^ 



+ ^ (2 - e'"J°^ 

90 #0 






(2^) 



T'Aoo(0,g) 



e ...,^) ^oo(go,g) ^^ (33) 



I'o 



which contains fewer terms than the Wilson loop of ( 13 ). 



* The crucial difference to potential models is that we do not in- 
vestigate the single point t = f}, but it is the fuU Euclidean time 
dependence of the gauge fixed correlator that is used to reveal 
the values of the potential. 



1. Coulomb gauge 

In Coulomb gauge, the HTL Euclidean gluon propaga- 
tor reads 



A^,.(go,9) 



cab 



Q' 



JQ{x-y) 



9fj.o9i^o 



Q2 + nT(Q) 



(?2 Q2 + ni(Q) 



(34) 



where the self energies TIl,t are the same as in covariant 
gauge (see Appendix Q . Inserting the propagator into 
the expression ( 33 ) for the Wilson line correlator gives 
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We now rewrite the HTL self-energies as spectral func- 
tions, use the formulas collected in Appendix [A| to per- 
form the sum over qo and carry out the angular integra- 
tions: 
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We find that the Coulomb gauge Wilson line correlator 
features a similar structure as the Wilson loop 



<'(r,r) = iy//j(r,r) + I^(l(r,r) 



II 



(37) 



the symmetric expression however being of much simpler 
form, depending only on the longitudinal HTL spectral 
function. At this point we can already anticipate that it 
is these terms present in both Wilson loop and Wilson 
line correlator, which contribute to the values of the po- 
tential. In particular the cusp singularity connected to 
the transverse spectral function identified in HI A3 c is 
absent from the above expression. 



2. Potential from the Wilson line correlator 

As in the case of the Wilson loop, a closed expression 
for the potential can be obtained using 

^"^(2)^•, .^ (38) 
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In the infinite time limit one can make use of 
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which leads us to the same result we encountered for the 
Wilson loop 
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with the imaginary part given by the integral expression 
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(41) 



From a practical standpoint this result is encouraging, 
as it tells us that (to leading order in HTL) the infor- 
mation content regarding the potential encoded in the 
Coulomb gauge Wilson line correlator is the same as the 
one found in the Wilson loop. If such a relation persisted 
into the non-perturbative realm, the absence of cusp di- 
vergences and with it the improved signal to noise ratio 
would make this an ideal observable to reconstruct the 
potential. 



3. Numerical evaluation 

As for the Wilson loop we wish to compare the Eu- 
clidean HTL correlator to actual values measured in 
quenched lattice QCD Monte Carlo simulations. While 
the symmetric term Wsym{T,r) is finite, the part linear 
in r still requires a regularization. We deploy the same 



momentum space cutoff as introduced in III A 2 and set 
its value to A = SttGcV in the following. 

The absence of divergences in the symmetric part of 
the correlator leads to a significantly different behavior 
along the imaginary times t. As can be seen in the top 
graph in Fig. |3J where we plot the HTL Wilson line corre- 
lator and the first five HTL Wilson loops as comparison. 
The large suppression at early times as well as the upward 
trend near r = /? are almost absent. Hence most of the 
datapoints actually carry information on the potential. 

Interestingly in the case of the lattice QCD Wilson 
line correlator, the upward trend is still visible between 
the last and second to last time step. However contrary 
to the leading order HTL result, where W^'^^{P,r) = 

W§'^^{/3, r) , the values of these two different correlators 
on the lattice do not agree at r = /?. 



4- Covariant gauge 

The Wilson line correlator can be calculated in the co- 
variant gauge as well. The result depends on the gauge 
parameter ^, and contains additional end point diver- 
gences [35j . These terms however do not contribute in 
the infinite time limit so that the obtained potential is 
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FIG. 3. (top) The Euclidean time Coulomb gauge HTL Wil- 
son line correlator VKi^"'^'^ (r, r) with momentum regulariza- 
tion A = Stt GeVevaluated at T = 2.33 x 270 MeVin steps 
of Ar = 0.066fm. (bottom) Quenched lattice QCD Wilson 
line correlator fixed to coulomb gauge from a lattice with 
a = 0.039fm and anisotropy ^ = 4 at T = 2.33Tc. Note 
that contrary to the HTL result the two correlators do not 
agree at t = /3 on the lattice. 



IV. SPECTRAL FUNCTIONS FROM HTL 
RESUMMED PERTURBATION THEORY 

A. Spectrum of the Wilson loop 

The spectral function can be directly calculated from 
the real-time correlator via a Fourier transform. 



Pui'T-.i^) = 
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2^ 



dt e«'^*es'^n '(**.'■) + 0{g^). (42) 



again the same as in the Wilson loop case ( 40 1 



We start by analytically investigating the low frequency 
behavior of this function, as it allows insight into the 
spectral structures that encode the physics of the heavy 
quark potential and will be used in its extraction in sec- 
tion [Vj To benchmark the MEM extraction of the spec- 



trum from Euclidean correlators it is however necessary 
to compare to the full spectrum, which we will determine 



from Eq. ( 42 1 numerically. 



1. Analytical estimate for the low energy part of the 
spectral function 



Starting form equation (21 1, we introduce the momen- 
tum cutoff A 
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where the argument of the second exponential function 
reads 
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For small freq uen cies, the main contribution to the spec- 
tral function (43) comes form small values of q" in the 



above integral. Expanding equation (45) around q^ = Q 
gives: 
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All terms with negative powers of q^ are retained in this 
expression, as they dominate the integral for late times. 
Note that the imaginary part of the potential appears as 
an overall factor in the above expression. Within this 
approximation, the remaining integrals are carried out 
analytically and we get: 
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(46) 



with d\j = '^"^^r'^''^' ^^'^ ^ denoting a not near specified 



_ |Im[y](r)| 
2T 

normalization constant. From this result, we see that 
the pole of the spectral function indeed resides at w = 
Re[y](r) and the width of the peak is closely related to 
the imaginary part of the potential. The result however is 
not a Lorentzian, but is precisely of the form (p| derived 
on general grounds in 15J. Note that the phase related 
to the skewing of the spectral peak is interestingly also 
given by the imaginary part of the potential 



Re[aoo]^5a = \MVKr)\/2T. 



(47) 
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FIG. 4. The HTL Wilson loop spectral function p^{r, uj) for 
different spatial separations Ar = 0.065fm. Note that the 
peak is extremely sharp but that its amplitude becomes very 
small at large r in comparison to the huge background induced 
mostly by the cusp divergences. 



2. Full spectral function 



We proceed to calculate the full spectral function by 
integrating numerically equation ( 42 ) . Applying the dis- 



crete Fourier transform to the real-time Wilson loop eval- 
uated on a set of Nt — 25000 points separated by a 



At 



50 GoV ' 



we obtain its values for a wide range of 



frequencies partly shown in Fig. [4] 



As expected from the minute values of fm[V ] at 
small separation distances, the peak one finds is ex- 
tremely sharp. However it also becomes clear that the 
amplitude of the peak is rapidly suppressed as r increases. 
At the same time non-potential contributions related to 
the divergent terms in the symmetric part of W^'^'{t,r) 
give rise to a huge background structure spanning a wide 
range of frequencies. 



Note that at oj « 18 GeV a step in the otherwise 
smooth spectral function is visible. This is a manifes- 
tation of the momentum cutoff we introduced to regular- 
ize the formally divergent terms. At the same time one 
can observe that the spectrum continues beyond these 
frequencies, which is a reminder that the cutoff was not 
imposed on the HTL gluon spectral functions. 



In section Iv] we will use the fitting function ^ to at- 
tempt an extraction of the heavy quark potential from 
the low frequency structures depicted in Fig. l4J 
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B. Spectrum of the Wilson line correlator in 
Coulomb gauge 

Analogously we can obtain the the spectral function 
related to the real time Wilson loop correlator 

1 



= ^ fdt e-^e^'^'^"'"^^ + 0{g^). (48) 
2tt J 

At leading order in the HTL resummed expansion, we 
again have: 
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The spectral function can then be calculated analytically 
close to its peak at small frequency, which yields 

A, ^ 1 |lm[T/](r)|cos(5|| - (Re[T/](r) - a.) sin,5|| 

^11^'^''^^ TT (Im[F](r))2 + (Re[F](r)-cj)2 

Surprisingly at leading order in the HTL approxima- 
tion we find that the skewing characterized by the quan- 
tity (5|| = ' ^^^r^^^'' is exactly the same as for the Wilson 
loop. Note that the same result can also be obtained in 
the covariant gauge. 



3, 




0.001 - 



FIG. 5. The spectral function of the HTL Wilson line correla- 
tor in Coulomb gauge Pu{r, to) for different spatial separations 
Ar — 0.065fni. While the peak position, width and skewing is 
exactly as in the Wilson loop case (Fig.p|, the absence of the 
cusp divergences leads to a significantly reduced background 
and a much higher amplitude at larger separation distances. 
Note that the plotting range is much smaller than in Fig. HI 




1. Full spectral function 

The full spectral functions for the HTL Wilson line 
correlator are plotted in Fig. [5J One immediately real- 
izes from a comparison with Fig. [4] that even though the 
peak position, width and skewing are equal to the Wil- 
son loop case, the Coulomb gauge spectral function looks 
quite different. The first major difference is that the am- 
plitude of the lowest lying peak depends much less on the 
separation distance r, the second is the virtual absence of 
the background terms populating a large frequency range 
in the Wilson loop case. Both facts are of course related, 
since their origin lies in the suppression of the Euclidean 
Wilson loop correlator induced in the presence of cusp 
divergences. 



V. THE POTENTIAL FROM PERTURBATIVE 
SPECTRAL FUNCTIONS 

Now that we arc in possession of the full HTL spec- 
tra obtained from both the Wilson loop and the Wilson 
line correlator in Coulomb gauge, we can test whether 



FIG. 6. Fits to the UV regularized (A = Stt) HTL 
spectral function pQ{r,uj) at r = 0.49fm (right) with 
a naive Lorentzian (L), a skewed Lorentzian (LS) and 
a skewed Lorentzian with additional polynomial terms 
(LSC0,LSC1,LSC2). Note that only the blue points (labeled 
"Fitted" ) are used for the fit and hence only these points enter 
the determination of the potential. 



the knowledge of the lowest lying spectral features alone 
suffices to reconstruct the values of the inter-quark po- 
tential in practice. To this end we fit the low lo re- 
gion of p^{r,uj) using the functional form ^ and com- 
pare the extracted values with the analytically calculated 
V^'^^{r). We show here the fitting of the Wilson loop 
spectrum only, since its application to yO||(r, w) gives ex- 
actly the same results (the potential and the skewing are 
the same) . In section |VI[ where the numerical recon- 
struction of the spectra from Euclidean time correlator 
data is concerned, the differences in e.g. the background 
contributions will however play a major role. 

In the following we do not constrain any of the pos- 
sible fitting parameters, i.e we allow e.g. Re[tToo] and 
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FIG. 7. (top) Real- and (bottom) imaginary part of the UV 
regularized (A — 57rGeV) HTL heavy quark potential (red, 
solid line) at T = 2.33Tc as well as the potential without cut- 
off (gray, dashed line). The various symbols denote the ex- 
tracted values from fits of the HTL Wilson loop spectra based 
on a Lorentzian (L), skewed Lorentzian (LS) and a skewed 
Lorentzian with background terms (LSCO, LSCl, LSC2). Note 
that the simple Lorentzian consistently overestimates the cor- 
rect values. The determination of the real part suffers only 
slightly from a worsening of the fit {LSC2 — >■ LS) but rough 
agreement is still visible. On the other hand a successful ex- 
traction of the imaginary part requires at least the presence 
of the first background term (LSCO), once the width of the 
spectral peak lies above 150MeV. 



Im[V] to be determined separately by the fit. To esti- 
mate in vifhich cases the use of improved fitting functions 
becomes necessary, we compare the results from a simple 
Lorentzian (L) (i.e. Re[croo] = and c^ = in eq.(|9|), 
a skewed Lorentzian (LS) (i.e. Ci — in eq.([9|) and 
the skewed Lorentzian with additional polynomial terms 
(LSCO, LSCl, LSC2) (i.e. c,>o = 0, c,>i = 0, c,>2 = 
ineq.([9|). 

We find that fitting with a simple Lorentzian (L) yields 
reasonable results only at very small separation distances 
r < 0.2fm, where the width of the peak itself does not 
exceed lOOMeV. In this distance region, the use of a 
skewed Lorentzian (LS) improves the fit significantly and 
actually reproduces the peak shape quite well. 

As shown in Fig. |6) at separations of r = 0.49fm 
the situation is already more involved, as the width of 
the peaks grows to around 150MeV and the shape de- 
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FIG. 8. Visualization of the connection 5u = 2t be- 
tween skewing parameter and imaginary part from the poten- 
tial peak in the HTL Wilson loop spectrum. For the skewing 
to be correctly determined, background terms up to quadratic 
order need to be taken into account. 



viates markedly from a naive Lorentzian. Adding the 
extra degree of freedom of skewing alone does not rem- 
edy the situation. Only after including the constant term 
(LSCO), arising from the early time variation of the func- 
tion $(r, t), we find that that the spectral shape is recon- 
structed in an acceptable manner. Including the addi- 
tional linear (ci) and quadratic (02) coefficients improves 
the overall agreement with the spectral shape, while the 
extracted values of the peak position and width are un- 
affected. This stability against including higher order 
background terms gives us confidence in the reliability of 
the fit. 

After scrutinizing the goodness of fit, we can turn our 
attention to the actual values of the potential obtained 
in this manner. In Fig. [7] the values of the real- and 
imaginary part of the UV regularized potential V (r) are 
shown in red (solid line) and the values obtained from the 
Wilson loop spectra fits are overlayed as discrete points. 
Wc find that the (LSC) fit successfully reproduces the 
real part of the potential at least up to the fourth digit. 
Note that the real part of the regularized potential shows 
an oscillating pattern absent in the UV complete V{R), 
which is retraced by the (LSC) fit. In lattice QCD where 
both a UV and IR cutoff are present, similar oscillations 
might arise. 

As expected from the fitting of the spectral shapes, the 
determination of the real and imaginary part succeeds 
even for the naive Lorentzian as long as the width is below 
lOOMeV. We find that the real part is less sensitive to the 
fitting function, only at larger distances, the Lorentzian 
overshoots the correct value, preventing us to observe the 
effect of Debye screening. 

The values of the imaginary part show a stronger de- 
pendence on the fitting function and the correct val- 
ues are only obtained after including the constant term 
(LSCO). In particular at large distances r > 0.3fm the 
simple Lorentzian and even the skewed Lorentzian over- 
estimate the values of the imaginary part. 



12 



According to the relation in Eq. ( 47 ) is should also be 



possible to extract the imaginary part of the potential 
through the skewing parameter Sa- While for a correct 
determination of the skewing, a precise fit of the peak 
shape is necessary, it is indeed possible to use the (LSC) 
fit to successfully relate skewing and imaginary part of 
the potential as shown Fig. [8j 

We conclude hence that the extraction of both the real 
and imaginary part from Wilson loop spectra succeeds 
if the improved fitting function eq. Q is deployecr] and 
therefore the knowledge of the shape of the lowest lying 
peak is sufficient to determine the potential. The results 
obtained with the (LSC) fit show a negligible deviation 
from the correct results and the deviation can be esti- 
mated by observing the variation of the fit results when 
introducing new fit parameters e.g. ci. We also find that 
fitting the lowest peak with a simple Lorentzian leads to 
an overestimation of both the real and imaginary part of 
the potential which contributed to the counter-intuitive 
results of Ref. [131 . 



VI. MEM ANALYSIS OF THE PERTURBATIVE 
EUCLIDEAN CORRELATOR 



While the extraction of both real and imaginary part 
of the potential from the lowest lying peak structure in 
p^{r, oj) has been shown to succeed in case of known HTL 
spectra in section Ivj we now wish to face the numerically 
challenging aspect of actually reconstructing these spec- 
tra from a set of Euclidean time data points. 

In the following we will deploy an MEM implementa- 
tion with extended search space [3B] (for technical back- 
ground see e.g. [57HiD] ) in an attempt to reconstruct 
from iV,- — 32 ideal imaginary time datapoints the most 
probable spectral function in the Bayesian sense. This 
number of available measurements along Euclidean time 
is representative for what we encounter in actual lattice 
QCD studies of correlation functions. By not adding ad- 
ditional noise and merely attaching artificial error bars 
to the correlators before feeding them to the MEM we 
deliberately choose the best case scenario in which any 
useful algorithm has to prevaij^ 



^ We checked that by including the next higher order in eq. (J9l, 
i.e. the term hnear in frequency with ci ^ 0, improves the fit 
at larger uj > 3.3GeV but does not change the extraction of 
the parameter values. If we go to higher temperatures, where 
the width becomes even larger or if we wish to fit the spectrum 
over a larger frequency interval around the peak, we will have to 
include higher terms of the c^'s. 

® One reasoning behind our choice is that e.g. through the ap- 
plication of the multilevel algorithm it is possible to measure 
datapoints with very high accuracy by sacrificing a number of 
available datapoints. 
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FIG. 9. MEM reconstructed spectra (dashed, blue) at r = 
0.066fm (top) and r = 0.264fm (bottom) based on iV^ = 32 
ideal Euclidean HTL Wilson loop data points at T = 2.33Tc- 
We discretize the frequency interval I^ = [—126, 189]GeV by 
N^ — 800 points and provide A^ext = Nr + 48 basis functions 
for the minimizer to reconstruct Pq^^{lu). The exact HTL re- 
sult at the corresponding distance is given as gray soHd curve. 
Note that even though the MEM recognizes the presence of 
the large background terms it fails to produce a smooth re- 
construction. Both position and shape of the lowest lying 
peak are rather poorly captured, which we attribute in part 
to the presence of the large background contribution. The 
limited number of available degrees of freedom do not suffice 
to capture both small and large {to > 5GeV) structures. 



2. Wilson loop 

To choose appropriate parameters for the MEM we 
first inspect the Euclidean data points in the top graph 
of Fig. [2J The strong suppression at small r as well 
as the rise at r ~ /3 tell us that structures at large 
positive and negative frequencies contribute to the full 
spectrum. Thus we decide to discretize w in an interval 
/„ = [-126,189]GeV by N^ = 800 points using arith- 
metic with a precision of 384 bits. 

The necessity for a large negative value of Wmin , indi- 
cated by the data, implies that the Nr basis functions in 
Bryan's search space do not contain enough variation to 
capture any peak at positive frequency. Hence we amend 
the search space by 48 additional basis functions of the 
full search space, whose oscillations cover the whole range 
of u!. The Levenberg-Marquardt (LM) algorithm is sub- 
sequently used to perform a search for the most probable 
spectral function within the confines of the above param- 
eters. 

Two of the resulting spectra are plotted in Fig. |9] We 
find that while the presence of the large background is 
acknowledged by the MEM through several peaks at fre- 
quencies 5GeV < a; < 50GeV, it is at the same time 
difficult to obtain a good reconstruction of the lowest ly- 
ing peak. At r = 0.066 at least its position is captured 
satisfactorily, the width of the structure remains an order 
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FIG. 10. Real part of the potential extracted from the MEM 
reconstructed HTL Wilson loop spectrum at T = 2.33Tc. 
We observe a consistent overestimation of the peak position, 
which persists even if higher background terms are included 
in the fitting function (e.g. LSC2). From the results of section 
[yjit is apparent that this failure originates in a deficiency of 
the underlying MEM reconstructed spectra. 
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FIG. 11. Imaginary part of the potential extracted from 
the MEM reconstructed HTL Wilson loop spectrum at T = 
2.33rc. We observe a consistent overestimation of the peak 
width of more than one order of magnitude. Note that in- 
cluding more fitting parameters worsens the estimation of the 
peak, since the underlying spectra do not actually resemble a 
skewed Lorentzian. 



of magnitude too large. 

Based on the MEM spectra we can proceed to fit the 
lowest peak using the fitting function Eq. ^ analogous to 
the spectra of section |V] The results are given in Fig. [TO] 
and Fig. [TT] The inadequacy of the spectral reconstruc- 
tion translates here into a consistent overestimation of 
the values for both real- and imaginary part of the po- 
tential. The shift in the peak position can be understood 
again from the presence of the large cusp divergence in- 
duced background, which together with the limited num- 
ber of basis functions, pulls the peak towards higher lo in 
the reconstruction. Similar to observations in previous 
lattice QCD based studies, both real and imaginary part 
are of the same order of magnitude. 

A few technical comments are in order. Our condition 
for finding the optimal spectrum in the Bayesian sense re- 
lies on a manual stopping criterion for the LM algorithm 
at relative improvements in the search of A = 10~^° to 
limit the necessary time for one run to the order of days. 
This however is not yet a true minimum since the val- 
ues meander around in tiny steps inside the search space 
without converging to a definite value within machine 
precision. This fact is sometimes refiected in non-smooth 
behavior of the a probability distribution. 

Increasing the number of basis functions improves the 
reconstruction slightly, i.e. the width decreases, but 
even with A'cxt > 200 we are not able to reproduce a 
Lorentzian peak shape including the characteristic tail 
structures. We furthermore only see marginal improve- 
ment in determining the lowest peak neither if the num- 
ber of datapoints is increased nor if the size of the ar- 
tificially attached errorbars is lowered. Since Bayesian 
inference is based on sound statistical reasoning with a 
well defined limit for infinitely many datapoints and ideal 
data these findings lead us to the conclusion that at this 



point it is not the properties of the supplied data but 
rather the implementation of the method that prevents 
us from a successful spectral reconstruction. 

It should however be noted that the difficulty in cor- 
rectly reconstructing the lowest peak is more difficult in 
the HTL case than what we expect to face on the lattice. 
Even though remnants of the momentum cutoff A are 
found in the HTL spectrum, the fact that higher frequen- 
cies contribute to integrals within the HTL gluon spec- 
tral functions, used at intermediate steps of determining 
^^^"'"'^(r, r), allows the background to stretch far beyond 
our choice of A = SttGcV. The presence of a sharp lattice 
cutoff would amputate such structures, the correspond- 
ing lattice correlator is less suppressed and thus the po- 
tential peak more easily reconstructed. 

We arrive at a sobering conclusion. Based on the Max- 
imum Entropy Method in its current form, even after in- 
cluding an extended search space, the reconstruction of 
the real and imaginary part from the Wilson loop is ex- 
tremely challenging. One of the reasons is the presence of 
the large background structures introduced by the cusp 
divergences (see Fig. |9| , which furthermore suppress the 
amplitude of the lowest lying peak. All attempts at a 
reconstruction of a sharp peak at small w are hampered 
since our limited reservoir of available degrees of freedom 
is depleted by structures not related to the physics of the 
potential. 



3. Wilson line correlator 

The reason to investigate alternative observables such 
as the Wilson line correlator in Coulomb gauge as basis 
for an MEM reconstruction is now evident. As we have 



seen in section III the absence of cusp divergences leads to 
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00 [GeV] 

FIG. 12. MEM reconstructed spectra (dashed, blue) at r = 
0.066fm (top) and r = 0.264fm (bottom) based on iV^ = 32 
ideal Euclidean HTL Wilson line correlator data points. We 
discretize the frequency interval lu, = [—63, 126] GeV by A'„ = 
2500 points and provide A^ext = ^t + 48 basis functions for 
the minimizer to reconstruct p?|'^'^(aj). The exact HTL result 
at the corresponding distance is given as gray solid curve. 
Note the absence of large background terms which hampered 
the reconstruction in the Wilson loop case. While the peak 
position is captured in a satisfying manner, the width of the 
reconstructed peaks is almost two orders of magnitude too 
large 



a dramatically reduced suppression along the Euclidean 
time axis. The rise at t = /3 observed in the Wilson loop 
is also virtually absent. This bodes well for an application 
of the standard MEM as the difficulties encountered in 
the previous subsection were directly connected with the 
divergences induced background contributions. 

We choose to discretize frequencies in an interval 1^ — 
[-63, 126]GeV by N^^ = 2500 points. The different choice 
of frequency range and N^^ compared to the Wilson loop 
case reflects our expectation that the available degrees 
of freedom suffice to reconstruct a much more narrow 
lowest lying peak. The necessity for accommodating a 
large background is gone. 



Fig. 12 shows two of the resulting reconstructed spec- 
tra which exhibit a much better agreement with the cor- 
rect HTL result than in the Wilson loop case. Note the 
factor five in the frequency axis compared Fig. [9] Car- 
rying out a fit with EqJ9] as in section |VJ allows us to 
estimate the values of the real- and imaginary part of 
the potential shown in Fig. [13] and Fig. [T4j 

We observe that at least for the real part a reason- 
able agreement with the correct Re[V^"'"'^] has been ob- 
tained, once the skewing is included (LS). What is strik- 
ing however is that the the values for different fit func- 
tions (LSCO, LSCl, LSC2) do not yet seem to asymptote 
for larger separation distances and thus begin to under- 
estimate the correct values. This behavior is connected 
to the fact that the shape of the reconstructed spectral 
peak does not resemble a skewed Lorentzian as can be 
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FIG. 13. Real part of the potential extracted from the MEM 
reconstructed HTL Wilson line correlator spectrum. After 
inclusion of skewing a reasonable agreement with the HTL 
potential is obtained. Note however that the values at larger 
r did not yet asymptote with respect to the inclusion of higher 
orders of background terms and tend to underestimate the 
correct values. The naive Lorentzian fit on the other hand 
leads to values that are too large. 
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FIG. 14. Imaginary part of the potential extracted from the 
MEM reconstructed HTL Wilson line correlator spectrum. 
Even though a much better resemblance of the reconstructed 
peaks with the exact HTL spectrum is obtained, we are still a 
factor five away from the actual values of the HTL imaginary 
part at this temperature. 



seen in Fig. [T2j 

The estimation of the width of the peak still fares worse 
compared to the peak position even though the disagree- 
ment has been roughly reduced by a factor three. The 
absence of the divergences in the Wilson line correlator 
and hence the absence of background terms already leads 
to a much more narrow reconstructed peak compared to 
the Wilson loop scenario, it is however still not possible 
to reach the actual width of the exact result. 

Since the reconstruction of the peak improved signif- 
icantly for the Wilson line case, we wish to inspect, 
whether the relation between skewing and peak width 
is already visible in the MEM spectra. Fig. [15] depicts 
the ratio between the fitting parameter (5)^'^'^ and Im[T^] 
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FIG. 15. Test of the relation (5|| — ™Jp ' between the the fitted 
skewing parameter and the imaginary part. We find signifi- 
cant deviations from unity with no clear tendency of improve- 
ment. This tells us that at this stage, the reconstructed spec- 
tral shapes still do not reliably reproduce a skewed Lorentzian. 



scaled by twice the temperature. The deviations from 
unity with no clear tendency of improvement tell us that 
at this stage, the reconstructed spectral shapes still do 
not reliably reproduce the skewed Lorentzian functional 
form actually encoded in the HTL Euclidean data. 

If the relation between skewing and imaginary part 
should turn out to hold beyond the leading order HTL 
approximation it would lend itself to checking the success 
of the MEM reconstruction. Note that in the Wilson loop 
case the extracted values of the skewing, besides being 
completely unstable between different fits did not show 
any correlation with the peak width. 

Despite the obvious technical shortcomings, which 
hamper the numerical determination of the potential, our 
findings are encouraging in that they show how the choice 
of underlying observable can improve the chances for a 
successful extraction of the potential. 

At least in the leading order HTL approximation the 
late real-time physics content is the same whether the 
Wilson loop or the Wilson line correlator in Coulomb 
gauge is concerned. The absence of divergences in the 
latter however permits the MEM reconstruction to lie 
much closer to the correct values. 

From the point of view of lattice QCD practitioners, 
the effects of e.g. the smearing procedure on spatial links 
in the Wilson loop are connected with modifying physics 
near the UV cutoff and might therefore also lead to an 
improved observable with respect to the potential extrac- 
tion. 



VII. CONCLUSIONS AND OUTLOOK 

The heavy quark correlator satisfies a Schrodinger 
equation at late real-times, parametrized by a complex 
potential. It has been proposed that this potential can 
in principle be extracted from the imaginary time Wil- 



son loop measured on the lattice [T31 [H] ■ The required 
steps involve an analytic continuation from Euclidean 
to Minkowski times, usually performed with the help of 
Bayesian inference and a fitting procedure that connects 
the position and shape of the lowest lying spectral peak 
to the real- and imaginary part of the potential respec- 
tively. 

Using HTL perturbation theory, we performed a sys- 
tematic check of the method with the aim to recover the 
well known potential of Ref. [5]. At leading order all 
quantities relevant to the extraction of the potential can 
be calculated explicitly. After the introduction of a mo- 
mentum cutoff for large spatial momenta, we were able 
to determine the full time dependence of both the Wil- 
son loop and Coulomb gauge fixed Wilson line correlator 
in the Euclidean as well as Minkowskian setting. The 
discrete Fourier transform is used to calculate the corre- 
sponding spectral functions. 

We find that the expectations of Ref. [TS] are fulfilled, 
as all spectra contain a well defined lowest lying peak 
of skewed Lorentzian form. A major difference between 
Wilson loop and Wilson line correlator in Coulomb gauge 
is the presence of cusp divergences in the former, which 
translates into large background structures engulfing the 
potential peak. Nevertheless the exact same potential is 
encoded in the two observables at this order in HTL. 

A surprising fact is that the skewing parameter, re- 
lated to the non-potential physics at early times, is the 
same for both observables and itself related to the imagi- 
nary part of the potential (see Eq.(40l). It would be very 
interesting to study both whether both of these proper- 
ties hold beyond the leading order HTL approximation 
or even on the lattice. 

Based on the HTL spectra we checked whether fitting 
the functional form derived in J15) succeeds and found 
that indeed a very accurate determination of the poten- 
tial is possible from exclusive knowledge of the lowest 
lying peak. We can hence replace late real-time limit in 
nh by fitting the low frequency realm of the spectrum. 
We find that the simple Lorentzian consistently overesti- 
mates the potential for intermediate and large values of 
the separation distance r, which offers a partial expla- 
nation for the large values observed in previous studies. 
By including skewing and higher order background terms 
from ([9]) however we obtain very stable and reliable esti- 
mates for Re[T^] and Im[F]. 

Our attempts of a standard MEM reconstruction of the 
spectra from ideal Euclidean time data-points revealed 
several challenges. We find that the large background, 
induced by the cusp divergences contained in the Wilson 
loop, makes a reconstruction very difficult. The broad 
structures not related to the potential absorb a large part 
of the limited number of degrees of freedom available to 
the MEM and prevent the potential peak to be captured 
in a satisfying manner. In addition we are not able to 
obtain a peak shape that resembles a skewed Lorentzian 
even after increasing the number of basis function for 
the MEM search space or adding additional datapoints. 
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This leads us to conclude that it is not the properties 
of the supplied datapoints but instead the technical im- 
plementation of the MEM itself that prevents us from a 
successful reconstruction. 

The situation is significantly improved in the case of 
the Wilson line correlator in Coulomb gauge, since the 
cusp divergences and hence the large background contri- 
butions in the spectrum are absent. Within the standard 
MEM it becomes possible to obtain a reasonable estimate 
of the real part of the potential, the imaginary part is still 
overestimated by at least a factor of five. It is still not 
possible to reconstruct the Lorentzian functional form 
encoded in the Euclidean correlator and thus the fits of 
([9]) to the both real and imaginary part become unstable 
if too many free parameters (e.g. q's) are included. The 
favorable UV structure of the Wilson line correlator in- 
vites speculation whether it can also help us to determine 
the non-perturbative potential from lattice QCD. 

At zero temperature, the potential is extracted from 
the Coulomb gauge Wilson lines anyway. In this case, it 
was shown in perturbation theory that at least to NNLO 
[IT] the full Wilson loop matches the Coulomb gauge Wil- 
son lines. The general picture is that if a physical poten- 
tial description exists at all in the large time limit, the 
details of the initial condition (how we close the Wilson 
loop) will not matteiN 

At non-zero temperature, if the Euclidean correlator is 
considered, there is no meaning to a large time limit as 
T < /3. However if we consider the analytically continued 
real time correlator or its spectrum, the infinite time limit 
might be considered similarly to the way it is performed 
at zero temperature. We showed here that the leading 
order HTL results agree one might expect, as in the zero 
temperature case, that the large time limit of the corre- 
lators do not to depend on the way we close them at the 
boundary. Should this assumption hold, then the poten- 
tial can be extracted from the Wilson line spectrum (or 
the smeared Wilson loop) as well, although not directly 
from the Euclidean data, for which the infinite time limit 
does not make sense. 



Appendix A: HTL spectral functions 



Inverting the relation (16), we have 
P(9°) = 



0, Aiqo ^ -*(q" + *0+)) - Aiqo ^ -liq'' - iO+)) 



2^ 

(Al) 

with A(go) ~ Qo + 'i^ + II(go, q). The spectral functions 
have a pole (gluon for the transverse, plasmon for the 
longitudinal) in the region q^ > q and a continuous part 
in the region q^ < q. They are antisymmetric in q^ and 
we restrict the formulas below to g° > 0. Explicitly, the 
electric spectral function reads 



PL{q°,q) 
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For the transverse spectral function, we have 

PT{q',q)'"^'7r6[fT{q")], (A3) 

/ \ 9>'3">0 

PT{q ,q) = 
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Appendix B: Convergence of the pole contributions 
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For the longitudinal spectral function, the solution to 
the delta function behave at large q'^ as 



ql 



l + 2cxp -2 



q^ +ml 



-D 



The full integral is then convergent because the factor 

«2 
1 



(ql)' 



is exponentially small. 

For the transverse spectral function, at large q, the 
pole is sitting at 



'' unless, of course, one chooses a divergent gauge, for instance 




D ' Of "flog"! 



(Bl) 
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The contribution from the poles is basically given by in- 
tegrating with respect to g" using the (5-function 

which gives in the transverse case a linear divergence. 

Appendix C: Cusp divergence 

We calculate here the divergence coming form the pole 
of the transverse HTL spectral function ( 28 ) in dimen- 



sional regularization and show that it matches the cusp 
divergences. 

For T = 0, /3 the function h{T,q^) vanishes and hence 
we have no contribution, note that there are no cusps in 
these degenerate cases. For t ^ 0, /3, one can decompose 
the factor 



nB{qT)H'r,qT) 



1 



r,rq°T _ pW-r)q°. 



where the second term lead to a convergent integral. 



Using equation (Bl), the remaining part of the inte- 



grand in (281 can also be decomposed as 

2 — ^™'-''^> — cos{qr) — qrSi{qr) 

TT . 

— -r + nnite 



1 



where the terms omitted are UV and IR finite. Using di- 
mensional regularization, the linear divergent term drops 
and only the r and r independent term survives. 



W, 



(2) 



pole,T 



{T,r) 



27r2e 



finite. 



(CI) 



The cusp multiplicative divergence arising from one angle 
7 reads |29l |30] 



-'cusp 



1 






(1 + (tt — 7) cot 7) 



In the Wilson loop case [28] with four angles 7 = 7r/2, we 
get as leading order contribution precisely the divergent 
part of (fell). 
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